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ABSTRACT 

Recently, several white dwarfs with very strong surface magnetic fields have been observed. In this paper we 
explore the possibility that such stars could have sufficiently strong internal fields to alter their structure. We 
obtain a revised white dwarf mass-radius relation in the presence of strong internal magnetic fields. We first 
derive the equation of state for a fully degenerate ideal electron gas in a magnetic field using an Euler-MacLaurin 
expansion. We use this to obtain the mass-radius relation for magnetic "^He, '^C, and ^^Fe white dwarfs of uniform 
composition. 

Subject headings: stars: white dwarfs - stars: magnetic fields - stars: interiors 



1. INTRODUCTION 

A number of w hite dwarfs with strong magnetic fi elds have 
been discover ed (Kemp et al. 197CI|; Putney 1995 , Schmidt 
& Smith 1995; [Reimers et al. 1996^ and exte nsively studied 
( Fordan 1992| ; Angel 1978; [Chanmugam 1992 and references 
therein). Surface magnetic fields ranging from about 10^ G to 
10' G have be en detected in about 50 (2%) of the « 2100 known 
white dwarfs (Jordan 1997 and references therein). As relics of 
stellar interiors, the study of the magnetic fields in and around 
degenerate stars should give important information on the role 
such fields play in star formation and stellar evolution. How- 
ever, the origin and evolution of stellar magnetic fields remains 
obscure. 

As early as Ginzburg (1964) and Woltjer (1964) it was pro- 
posed that the magnetic flux ($5 ^ BR^) of a star is conserved 
during its evolution and subsequent collapse to form a remnant 
white dwarf or neutron star. A main sequence star with ra- 
dius on the order of R ^ 10^' cm and surface magnetic field 
B ^ 10- 10"* G [magnetic A-type stars ha ve typical surface 
fields < 10** G ( [Shapiro & Teukolsky 198^ )] would thus col- 
lapse to form a white dwarf with R ^ lO'^ cm and B ^ 10^-10^ 
G, or a neutron stai- with 7? - 10*' cm and B - 10" - 10'** G. In- 
deed, shortly after their discovery (Hewi sh et al. 19 68) pulsars 



were identified as rotating neutron stars (Gold 1968) with mag- 
netic fields B ^ 10"-10'"'G consistent with magnetic field am- 
plification by flux conservation. In addition, neutron stars with 
surface magnetic fields exceeding lO''* G [so called magnetars] 
have been recently suggested as the so urce of soft gamma-ray 
repe aters (0uncan & Thompson 1992|; T hompson & Duncan 
1995). 

Moreover, the surface magneti c field of a star d oes not nec- 
essarily reflect the internal field (Ruderman 1980). For exam- 



ple, the toroidal fields below the surface of the Sun are at least 
on the order of ~ 10^ to ~ 10** times stronger than the aver- 
age surface dipole field strength of ~ 1 G (G alloway, Proc- 
tor, & Weiss 1977). Furthermore, at the region of the con- 
vective zone, the strength of small scale magnetic fields could 
reac h a value as high as 7 x 10"* G (C hauhan, Pandey & Pandey 
1999, pulido 1998). This would correspond to an interior field 



strength on the order of ^ 10 to ^ 10 G in a white dwarf, 
or ^ 10'^ to ^ 10'^ G in a neutron star [Condensed objects of 
size R and mass M have an upper limit to their field strengths 



of B < M/?^(87rG)'''^. For neutron stars with /? w 10 km and 
MkjMq, the limit is B <^ lO''^ G (pLerche & Schramm 1977|).] 



Indeed, the existence of white dwarfs with interior magnetic 
fields as strong as ^ 4 x 10'^ G is not ruled out with the present 
unce rtainties in the mass-radius relation ( Shapiro & Teukolsky 
1983). The present high upper limit on the strength of internal 
fields in white dwarfs is obtained by simply setting the mag- 
netic pressure equal to the internal pressure of the star. How- 
ever, white dwarfs wit h internal field s at or around this strength 
could be constrained (Mestel 1965) by a perceptibly different 
mass-radius relation. 

Although white dwarfs in binaries with well determined 
masses do not appear to have surface magnetic fields larger than 
^ 10^ G, interna l fields of ord er lO'^ G could be well hidden be- 
low the surface (Angel 1978). Newly discovered magnetic de- 



generate stars, especially those with surface field strengths near 
the range of B ~ lO' G, always show strong circularly and/or 
linea rly polarized spectral energy distributions (S chmidt et al. 
1999). Moreover, these s tars reveal unique spectral features 



(Engelhardt & Bues 1994) due to quasi-Landau resonances in 
extremely high magnetic fields of > lO' G. 

In this work, we explicitly compute the mass-radius rela- 
tion of white dwarfs with internal magnetic fields. Previously, 
Ostriker & Hartwick (1968) have estimated effects of interior 
magnetic fields by considering a correction in terms of the ratio 
of magnetic to gravitational energy. They showed that a rel- 
atively small ratio of magnetic to gravitational energy would 
be sufficient to explain an observational discrepancy in the 
classical mass-radius relation for Sirius B. However, if white 
dwarfs could indeed have central magnetic fields as strong as 
4.4 X 10" -4.4 X 10'^ G, the revised mass-radius relation must 
be explicitly determined by taking the magnetic field into ac- 
count in the equation of state. The present work thus expands 
upon that earlier study by explicitly computing the equation of 
state for a completely degenerate, noninteracting electron gas 
in a magnetic field. This equation of state is then applied to the 
Tolmam-Oppenheimer-Volkoff (TOV) equation of stellar hy- 
drostatic equilibrium. 

The equation of state in a magnetic field should reduce to 
a normal equation of state in the absence of a magnetic field. 
Therefore, we use an Euler -MacLaurin expansion (K ernan, 
Starkman, & Vachaspati 1996) of the thermodynamic variables 
to recover the weak field limit. In integrating the TOV equation 
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we simply follow the procedure of Hamada & Salpeter (1961) 
for degenerate matter of uniform composition of ^He, '^C, or 
^^Fe. Although '^He and '^C white dwarfs are expected to have 
a similar (though not identical) mass-radius relation, we explic- 
itly consider each for completeness. 

2. EQUATION OF STATE FOR AN ELECTRON GAS IN A MAGNETIC 

FIELD 

The properties of an electro n in an external magneti c field 
have been studied ext ensively ([Landau & Lifshitz 1938|; John - 
son & Lippmann 1949; [Canuto & Chiu 1968| ; [Schwinger 1988[ ). 
In brief, the energy states of an electron in a magnetic field are 
quantized and its properties are modified accordingly. In order 
to investigate these effects, we must first solve the Dirac equa- 
tion in an external, static, and homogeneous magnetic field. We 
make the convenient choice of gauge for the vector potential 
in which a uniform magnetic field B lies along the z-axis. We 
then obtain the el ectron wavefunctions and ene rgy dispersion in 
a magnetic field (Johnson & Lippmann 1949). The dispersion 
relation for an electron propagating through a magnetic field is 



E„ = [p^c^ + n?c'^ + 2nceBn\ 



1/2 



(1) 



where n = j+ j 



+ s,,(n = 0,l, ...), j is the principal quantum 
number of the Landau level, 5, = ±1/2 is the electron spin, e is 
the electron charge, c is the speed of light, h is Planck's con- 
stant, p = pz is the electron momentum along the z-axis, and 
m is the rest mass of the electron. In Eq. (1) we ignore the 
anomalous magnetic moment for an electron. 

The main modification of an electron in a magnetic field 
comes from the available density of states for the electrons 
(Landau & Lifshitz 1938). The electron state density in the 



absence of a magnetic field. 



is replaced with 



d^p 

(2^ 



9 PR 

in a magnetic field. The symbol 5„q is the Kronecker delta de- 
fined by (5„o = 1 for n = 0, or (5,!0 = for « ^ 0. This modification 
affects the thermodynamic variables for the electron gas. 

An isolated white dwarf ultimately cools to zero temperature 
(fully degenerate) and the degeneracy pressure supports these 
stars against further gravitational collapse. For the most part, 
this pressure can be described as an ideal (noninteracting) elec- 
tron gas plus small corrections. 

Let us consider a gas of electrons at zero temperature in a 
magnetic field. From Eq. (1) we can define the Fermi energy 
Ep for an arbitrary Landau level n as 



E} = m^c'^ +pj.c^ + ITiceBn . 



(2) 



Here pf denotes the Fermi momentum. The number density of 
electrons in a magnetic field is then given by 



"^ = '(2^(t)'^^^^'"^' 



(3) 



where 



«/ 



C(e/,«) = ^[2-<5„o] Je)-{l+2^n) 



In the above, e/ = Ef/mc^ and 7 = B/Bc where Be = 
np-c^ /\e\Ti = 4.414 X lO'-' G is the critical magnetic field at 
which quantized cyclotron states begin to exist. The maximum 
Landau level nj for a given Fermi energy ey and magnetic field 
strength 7 is given by 
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The pressure of an ideal electron gas in a magnetic field is then 

3 



7 



o / mc \ ^ 



(4) 



where 



1 "' 
$(£/,«) = -Y^[2-S„o] 



n=0 



e/./e2-(l+27«) 



-(l+27n)ln 
Similarly, the energy density is 



'e/+^e2-(l+27«)^ 



£{ef,n) = 2-^mc (^— j x(e/:«), 



(5) 



where 



1 "' 
X(^f,n) = -^[2-(5„o] 



efJe}-(l+2jn) 



+(l+27n)ln 



'ef+^/ej-(l+2jnr 



^/l+2"fn 
From these, we obtain the energy per electron 

„2 X(e/,») 



Ee{ef,n) = mc 



C(<^f,n)' 



(6) 



In order to recover the usual equation of state in the absence 
of a magnetic field, we utilize an Euler-MacLaurin expansion of 
Eqs. (3) - (6) in the weak field limit. Then, the number density 
is given by 



where. 



"^"3^(t)'^^"^' 



ax)=x'+-f^-+o(j'')+- 



(7) 



and X = pf/mc is the relativity parameter. Note that the elec- 
tron number density increases as the magnetic field increases 
for a given x. 

The pressure becomes 



1 -) / mc \ 3 

r 



24.2-"^^nTJ *^^)' 



(8) 



where 



n=0 



<^ix) = <^o(x) + Y^Bix) + 0{J*) + • • • , 
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%(x) = xVx^+l(2x^-3)+3ln(x+\/x^+l), 



Jx^+1 , 

$g(x) = + 21n(x + Vjc2+l) 



X 

1 + - 



ix+Vx^ + 1) 



Note also that for a physically reasonable value of x the pressure 
always increases as 7 increases. 

The energy density can also be written 



where 



2 / "^c \ ■* 



X(x) = Xo(x) + 7^Xb W + 0(7"^) + • • • , 



(9) 



Xoix) = — -r 



x(2x^+ l)Vx2+T-ln(x+ V?+T) 



such as the Thomas-Fermi, exchange, and correlation correc- 
tions give only a very small change in the mass-radius relation 
of white dwarfs. Actually the Thomas-Fermi correction in a 
strong magnetic fiel d is important at low density (see Rogn- 
valdsson et al. 1993 and references therein), but as a whole it 
gives only a minor effect on the mass-radius relation. Hence, 
we ignore these minor effects for the present work. 

Magnetic fields should not alter the spherical symmetry of 
the constituent atoms. Indeed, it is, perhaps, a remarkable fact 
that magnetic fields do not destroy even the approximate spher- 
ical symmetry of he avy atoms within the relevant range ( Rogn- 
valdsson et al. 1993). Thus, we can use the ordinary electro- 
static energy (Coulomb energy) per charge. 



5 re 



(12) 



where Ry = ^a^nieC^ is the Rydberg energy and from Eq. (7), 



Xb(x) 
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Vx^TT 
1 + + - 



x(x+Vx^+l) 



-21n(x+Vjc2+l) 



Finally, the energy per electron is given by 



Eeix) ■ 



X(x) 



(10) 



Here we can see explicitly that as 7 goes to zero, Eqs. (7) - (10) 
recover exactly the usual equation of state in the absence of a 
magnetic field. 

3. MASS-RADIUS RELATION OF MAGNETIC WHITE DWARFS 

The mass-radius relation of white dwarfs was first deter- 
mined by Chandrasekhar (1939). Later Hamada & Salpeter 
(1961) obtained numerical models for various core composi- 
tions by considering a fully degenerate configuration at zero 
temperature. The theoretical relationship between the mass and 
radius of a white dwarf i s important for the interpretati on of 



observational results (see Koester & Chanmugam 1990 for a 
review). There are severa recent studies and observations on 
the mass-radius relation of non-magnetic white dwarfs (W ood 
1990; [Vauclair, Schmidt, & Koester 1997|; Vennes, Fontaine, & 
Brassard 1995; Provencal et al. 1998^ 

In order to obtain the mass-radius relation for magnetic white 
dwarfs, we use Eqs. (7) - (10) for 7 < 1 (B < 4.4 x 10'^ G) and 
carry out stellar integrations for a uniform composition of ^He, 
' ^C, and ^''Fe as an illustra ti ve model. In this un iform model 
( [Hamada & Salpeter 196^; pushiki et al. 1992"|; R ognvalds- 
son et al. 1993) the total energy E of the plasma consists of 
a nearly uniform distribution of degenerate electrons with em- 
bedded ions, 

E = Ee + Ec, (11) 

where the first term is the energy of a uniform gas of free elec- 
trons and Ec corrects for the classical Coulomb energy. Al- 
though the noninteracting electron gas accounts for the dom- 
inant contribution to the equation of state at high density, the 
classical Coulomb correction is significant. Other corrections. 



37r 



87CW 



1/3 



The corresponding pressure is 



, /OTc\3 9 Z^^a^ 1 

Pc = -mc^{ — J. 

V fi / 47r 10 rf 



(13) 



In this model, the total pressure is then given by 



P = Pe + Pc. 



(14) 



Thus, since Pc is negative, the eq uation of state would lead to 
negative pressures at low density ( Salpeter 1961 ). 

In integrating the equations of hydrostatic equilibrium, we 
follow the classical procedure of Hamada & Salpeter (1961). 
Figures 1, 3, and 5 show the mass-radius relation of "^He, '^C, 
and ^^Fe white dwarfs for a given magnetic field strength. Fig- 
ures 2, 4, and 6 show the relation between mass and cen- 
tral density of these white dwarfs for a given magnetic field 
strength. Here it can be seen that our results approach the 
Hamada & Salpeter (1961) results as the magnetic field strength 
decreases. For high central field strengths 7 ~ 0.01 - 1 [B ~ 
4.4 X (10" - lO'"*) G], both the mass and radius of magnetic 
white dwarfs increase compared to non-magnetic white dwarfs 
of the same central density. For instance, for 7 « 0.8 carbon 
white dwarfs, the radius increases by about 30% forM w IMq. 
Similarly for R w O.O17?0 the mass also increases by about 
25%. These results are approximately consistent with Ostriker 
& Hartwick (1968), i.e. the radius increases while the central 
density p^ decreases as the magnetic field increases for fixed 
M. As expected, forB < 10'*' G, internal magnetic fields do not 
affect the white-dwarf mass-radius relation. 
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Fig. 1 Relation between the mass M and radius R of a ''He magnetic 
white dwarf for the indicated magnetic-field strengths. The solid line 
denotes the Hamada & Salpeter model for non-magnetic white dwarfs 
('y = 0). The dashed lines are magnetic white dwarfs. 
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Fig. 2 Relation between the central density pc (in g/cm^) and mass M 
for ''He magnetic white dwarfs. The solid line denotes the Hamada- 
Salpeter model for non-magnetic white dwarfs (7 = 0). The dashed 
lines are for magnetic white dwarfs. 



A striking feature of these results is that white dwarfs with 
strong interior magnetic fields should be massive. This is sim- 
ply because a star becomes unbound if the magnetic plus mat- 
ter pressure force exceeds the gravitational force. For example, 
there is no stable solution with M <Q.5Mq for 7 = 0.3 carbon 
white dwarfs. This is consistent with recent observations that 



on average magnetic white dwarfs have a higher mass than typ- 
ical non-magnetic white dwarfs (see Table 1). Note, however, 
that for a field strength of B < 1 x 10'^ G, magnetic fields only 
give a relatively small change on the mass-radius relation, that 
is, we can not distinguish between magnetic and non-magnetic 
white dwarfs for 7 < 0.2. 

4. DISCUSSION 

In this work, we have calculated the equation of state for an 
electron gas in a magnetic field at zero temperature. For sim- 
plicity, we have assumed a uniform composition to obtain the 
relation between mass and radius for magnetic white dwarfs. 
For high internal magnetic fields B ~ 4.4 x (10^' - 10'^) G, 
(7 ~ 0.01 - 1), the mass-radius relation is modified. Our results 
not only confirm the Ostriker & Hartwick (1968) result of in- 
creasing radius and decreasing central density pc with increas- 
ing field, but also are consistent with the suggestion (Liebert 
1988) that observed magnetic white dwarfs have masses which 
are on average larger than non-magnetic white dwarfs, imply- 
ing more massive and younger progenitors. 

The question remains, however, as to whether it is reason- 
able to consider such high internal field strengths for magnetic 
white dwarfs. First, one must assume that the magnetic fields 
are well hidden beneath the surface, while the surface fields are 
several orders of magnitude less. Second, assuming that flux is 
conserved during the collapse to a white dwarf, the progenitor 
of the white dwarfs must have had sufficiently large fields to 
produce the required white dwarf internal field strengths. Flux 
conservation implies that the central field strength of the pro- 
genitors is of order ^ 10*^ G, assuming R k, \Rq. This is rea- 
sonable. During star formation, the collapse of a typical inter- 
stellar cloud with radius ^0.1 pc, mass ~ \Mq, and protostel- 
lar magnetic field of magnitude ~ 3 x 10~* G would result in a 
field strength of ^ 3 x 10^ G in a solar type star formed from 



this material ( Spitzer 1978 ). Although there is no evidence for 
main sequence stars with s uch field strengths, they are n ot ruled 
out by observations either ( Shapiro & Teukolsky 1983 ). 



5. COMPARISON WITH OBSERVATIONS 

Figure 7 shows a comparison between our calculations and 
white dwarfs with known masses and radii from the HIPPAR- 



COS survey ( [Vauclair, Schmidt, & Koester 1997| ; Provencal et 
al. 1998). Plotted error bars are the quoted ±lcr observational 
uncertainties. Strong hidden interior magnetic fields would be 
expected to manifest themselves by a preponderance of stars 
with large masses and radii. Most of the data, however, are 
within 2a of the non-magnetic theoret ical curves. A puzzling 
feature however ( Provencal et al. 1998) is that some of the best 



determined data points (e.g., EG 50 and Procyon B) can not be 
fit without postulating an iron composition, something which 
seems unlikely from a stellar evolution standpoint. Further- 
more, although the evidence is not compelling, there are at least 
two well d etermined stars [i.e., GD 140 and Sirius B (P rovencal 
et al. 1998)] as well as some field stars with M > Q.6Mq and 
R > Q.Q12Rq which may be better fit if strong internal mag- 
netic fields are assumed. However, Sirius B and some of the 
field stars ca n also be accommodated by atmospheric models 
(IWood 19901). 
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Fig. 3 Same as Fig. 1, but for '^C. 



Fig. 5 Same as Fig. 1, but for ^^Fe. 
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Fig. 4 Same as Fig. 2, but for '^C. 



Fig. 6 Same as Fig. 2, but for ^^Fe. 



Mass-Radius Relation for Magnetic White Dwarfs 




0.6 0.8 1.0 

M/Mq 

Fig. 6 A comparison of white dwarfs with known 
from the HIP^RCOSsurvevatthelaleveWyauclair, 
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1, 3, and 5. The upper curves are for He and C wi 
lower curves are for ^^Fe white dwarfs. 



Some typical magnetic white dwarfs which have reported 
masses M and surface magnetic fields Bs are summarized in 
Table 1 . This table shows that magnetic white c warfs are typ- 
ically more massive than non-magnetic white d\ /arfs (on aver- 
age M - O.6M0). Among them RE J03 17-858 {B, > 10*^ G) 
is approaching the Chandrasekhar mass limit (Barstow et al. 
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1995). 

Nc t listed in Table 1 , however, are the two known magnetic 
white dwarfs with the strongest surface magnetic field, GD 229 
(B, > 10'' G)andPG1031^-234(B,«(0.5-l)xlO^G).Un- 
fortu lately, there are no reported masses or radii for those mag- 
netic white dwarfs. If their masses and radii could be measured, 
then )ur magnetic field model might be tested. 

Pe "haps the most interesting object in Table 1 is LB 1 1 146 
(PG )945h-245). This is an unresolved binary sy stem consist- 
of two degenerate stars ( Liebert et al. 1993 ): one compo- 
s a normal DA white dwarf (M = 0.91 ± 0.07 Mq) with no 
detecltable magnetic field; the other has a strong magnetic field 
3 X 10*^ G) and a mass in the range 0.76 < M/Mq < 1 .0. 
, assuming for illustration a strong interior field of 7 ~ 0.5, 
that this is a typical equal mass system with M = 0.9 Mq 
then the radius of the magnetic component would be 
■ by about 10% than the normal star Alternatively, if they 
similar radii of R = 0.01 Rq, the magnetic star would have 
heavier mass by about 12% than the normal one for 7 ~ 0.5. 
Re garding the radius of magnetic white dwarfs, Greenstein 
01 :e (1982) have reported radii of 7? - 0.00667^0 for Grw 
B247 (B, « 3.2 x 10** G) and R - OMRq for Eeige 7 
^ 3.5 X 10^ G) from the interpretation of their spectra. 
Assuning an interior field which is about 10^ times stronger 
the surface field, then we would deduce 7 w 0.8 and 
1.1 - I.2M0 for Grw H-70°8247, consistent with the ex- 
pects tion of high masses for magnetic white dwarfs. Similarly, 
we w ould deduce 7 « 0.075 and M ^ O.8M0 for Eeige 7. This 
star lowever would be nearly indistinguishable from a non- 
magr etic white dwarf. 
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Mass and surface magnetic field strength in some typical magnetic white dwarfs 
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